Highly anisotropic resistivities in the double-exchange model for strained manganites 
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The highly anisotropic resistivities in strained manganites are theoretically studied using the two- 
orbital double-exchange model. At the nanoscale, the anisotropic double-exchange and Jahn- Teller 
distortions are found to be responsible for the robust anisotropic resistivities observed here via Monte 
Carlo simulations. An unbalanced in the population of orbitals caused by strain is responsible for 
these effects. In contrast, the anisotropic superexchange is found to be irrelevant to explain our 
results. Our model study suggests that highly anisotropic resistivities could be present in a wide 
range of strained manganites, even without (sub)micrometer-scale phase separation. In addition, 
our calculations also confirm the formation of anisotropic clusters in phase-separated manganites, 
which magnifies the anisotropic resistivities. 

PACS numbers: 75.47.Lx, 71.70.Ej, 75.30.Gw 



I. INTRODUCTION 

Strongly correlated electronic materials, which are 
well known for the presence of complex phase compe- 
titions involving the spin, charge, and orbital degrees of 
freedom,^ are promising candidates to be used in new 
multifunctional devicesi^ Typically, in materials such as 
manganites with the colossal magnetoresistance (CMR), 
there are several phases with free energies that are quite 
close to one another but their individual physical prop- 
erties can be rather different.'^ Therefore, colossal re- 
sponses to external perturbations, including the CMR*^ 
and colossal electroresistance (CER),^ can and do occur 
in some manganites. During the past decade, theoreti- 
cal studies on manganites have addressed many of these 
colossal responses, such as CMR,^"^ CER,^"'^^ surface 
reconstructions r^"— and disorder effects J^"— 

In addition, the effects of strain on the properties 
of manganites and other complex oxides is attracting 
increasing attention due to the rapidly expanding re- 
search interests in complex oxides heterostructures i^^i^^ 
In fact, phase transitions driven by strains have been dis- 
cussed in manganite thin films for several years.— ^— The 
physical mechanism of these phase transition is mostly 
orbital-order-mediated)^— For example, according to 
density functional theory (DFT) calculations, the ground 
states of LaMnOa/SrMnOs superlattices can be tuned 
between A-type antiferromagnetic, ferromagnetic, and C- 
type antiferromagnetic phases when the ratio c/ a is in 
the range 0.96 ^ 1.04, where c (a) is the out-of-plane 
(in-plane) lattice constant.-— Even for LaMnOa itself, 
the ground state may become ferromagnetic (FM) if the 
|3x^ — r^ >/|3y^ — r^ > type orbital order is fully sup- 
pressed in the cubic lattice, according to both the DFT 
and model calculations j^^*^ 



Very recently. Ward et al. have ob- 

served high anisotropic resistivities in strained 
Lag/s-^rPr^Cag/sMnOa (LPCMO) thin films.i^ LPCMO 
is a prototype phase-separated material4i The coex- 
istence of FM and charge-ordered-insulating (COI) 
clusters at the (sub) micrometer-scale can seriously affect 
the electric transport properties, especially the metal- 
insulator transition (MIT). The electric conductance 
in the phase-separated LPCMO is dominated by the 
percolation mechanism.— i^ For example, giant discrete 
steps in the MIT and a reemergent MIT occur in an 
artificially created microstructure of LPCMO when the 
size confinements in two directions become comparable 
to the phase-separated cluster sizes.— i^ Therefore, 
Ward et al. proposed that the anisotropic percola- 
tion might be responsible for the highly anisotropic 
resistivities in strained LPCMO. ^° Also, our previous 
simulation of CER predicted anisotropic resistivities 
due to the electric-field-driven anisotropic percolation in 
phase-separated manganitesJ^ 

Then, two interesting questions arise: (1) how does 
strain drive the anisotropic percolation in the LPCMO 
films? And, more importantly, (2) can the large 
anisotropics occur in more standard CMR materials with 
nanometer-scale phase competition or even with bicrit- 
ical clean-limit phase diagrams? Therefore, to setup a 
study to be used as a reference for future research it is 
interesting to investigate theoretically with model Hamil- 
tonians the magnitude of the anisotropy in transport 
induced by strain in cases where phase competition is 
present, but also where phase separation is not. In other 
words, it is important to study regimes where in the clean 
limit (no quenched disorder) a first-order transition sep- 
arates the two competing states, typically a metal and 
an insulator, inducing a CMR effect in a narrow range of 



parameters, but where phase separation it not present. 
This calculation will allow us to disentangle the effects 
of mere strain on a clean limit model in the regime of 
phase competition from the effects of strain on a truly 
phase separated state. More basically, these investiga- 
tions are important to move beyond the micrometer-scale 
to find the microscopic origin of anisotropic resistivities 
in generic strained manganites. 



II. MODELS AND TECHNIQUES 

In this paper, the two-orbital double-exchange (DE) 
model will be employed to study the anisotropic re- 
sistivities in strained manganites. In the past decade, 
the DE model has been extensively studied and it 
proved to be a quite reasonable model to describe per- 
ovskite manganitesi^ In Ward et al.'s experiments, the 
anisotropic strain field splits the in-plane lattice con- 
stants along the [100] and [010] axes in the pseudocubic 
convention (or the [101] and [101] axes in the orthorhom- 
bic Pnma convention). Thus, a modified model has to be 
developed to reflect the features of this strained lattice, 
since most previous model studies were done on cubic or 
square lattices. 

As a well-accepted approximation for manganite mod- 
els, an infinite Hund coupling is here adopted. With this 
useful simplification, the DE model Hamiltonian reads: 

a/3 
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In the above model Hamiltonian, the first term is 
the standard DE interaction, a and /? denote the two 
Mn Cg-orbitals a (|x^ - y^ >) and h (|3z^ - r^ >). 
Cia (Cict) annihilates (creates) an Cg electron at orbital 
a of site i, with its spin parallel to the localized i2g 
spin Si. The nearest- neighbor (NN) hopping direction 
is denoted by r. The Berry phase J7y generated by 
the infinite Hund coupling equals cos{9i/2) cos{6j/2) + 
sm{6i/2) sm{6j /2) exp[—i{(t)i — 4>j)], where 9 and are 
the polar and azimuthal angles of the i2g spins, respec- 
tively. In strained manganites, an elongated lattice con- 
stant gives rise to more straight Mn-O-Mn bonds, thus 
enhancing the FM DE interaction. To mimic this effect, 
the in-plane DE hopping amplitudes t^g have to be set 
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In the rest of the manuscript, ip is taken as the energy 
unit tg and At = t^/t^ — 1 is defined to characterize the 
degree of anisotropy of the DE interaction. 

The second term of the model Hamiltonian is the an- 
tiferromagnetic (AFM) superexchange (SE) interaction 
between NN t2g spins. The SE coefficient Jaf could also 
become anisotropic in the strained lattices, which is here 
characterized by Aj = J^f/'^af ~ ^■ 

The third term of the model stands for the electron- 
lattice coupling. A is a dimensionless coefficient and n^ 
is the eg electronic density at site i. Qs are phonons, in- 
cluding the Jahn- Teller (JT) modes {Q2 and Q3) and 
the breathing mode (Qi): Qi = {Sx + Sy + 6z)/V^, 
Q2 = {5x-5y)/^/2, andQa = (-J^, -(5j^ + 2(52 )/V6, where 
5 stands for the length change of the oxygen coordinates 
in the Mn-O-Mn bonds along the axes directions, r is the 
orbital pseudospin operator, namely Tx = c\ch + clca and 
Tz = c\ca — c\ch. The last term is the lattice elastic en- 
ergy. Note that the model used here induces cooperative 
distortions of the oxygen positions. 

The above model Hamiltonian is numerically solved via 
the Monte Carlo (MC) simulation on a two-dimensional 
8x8 lattice. The reason for this restriction to a two- 
dimensional geometry is simply practical: simulations in 
three-dimensional lattices are very demanding computa- 
tionally. Thus, here 5z is set to zero and our effort will 
only focus on the in-plane anisotropy. Using standard 
periodic boundary conditions (PBCs), < Qi >, < Q2 >, 
and < Qa > (if <> stands for averages over the whole 
lattice) equals to zero. However, to simulate the strain 
effect in the JT distortion, anisotropic PBCs (aPBCs) 
should be introduced to the lattice. In the aPBCs for 
2D lattices, < Q2 > is set as a constant which can be 
nonzero, while < Qi > and < Qz > remain zero. To 
characterize this anisotropic JT distortion, the quantity 
Aq is defined as - < (32 > 7(2^3). 

In Ward et aVs experiments, the difference between 
the in-plane lattice constants is small (~ 0.2 — 0.3%).— 
Correspondingly, the anisotropics of interactions should 
be weak, implying that At, Aj, and Aq must be small 
quantities in our study. 

In our MC simulations, the average Cg density < n > 
is chosen as 0.75. As discussed in previous literature, to 
obtain the MIT and CMR effects, the parameters (Jaf, 
A) should be chosen to be near the phase boundaries be- 
tween FM and AFM COI phases. ^'^ This fine tuning of 
couplings could be avoided by introducing quenched dis- 
order, but our study will be conducted in the clean limit 
to setup a benchmark to decide on the origin of strain 
induced transport anisotropics that are investigated ex- 
perimentally. According to the phase diagram of the 
two-orbital DE model for < n >= 0.75,^^ the parame- 
ters Jaf — 0.09 and A — 1.2 are suitable and they are 
here adopted as the default ones in our simulation, un- 
less other parameters are explicitly used. In fact, other 
sets of parameters near the default ones have also been 
partially tested and no qualitative differences have been 
found. Thus, this choice of parameters do not alter the 
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FIG. 1. (Color online) MC simulated resistivities (triangles) 
and magnetization (dots) for a square 8x8 isotropic lattice 
[At = 0, Aj = 0, and Aq = 0), as a function of T. 



general validation of our results and conclusions, at least 
qualitatively. In the MC simulation, the first 10'* MC 
steps are used to reach thermal equilibrium and another 
2 X 10^ MC steps are used for measurements. 

The dc conductances, which are calculated using the 
Kubo formula, are in units of e^//i, where e is the el- 
ementary charge and h is the Planck's constant^S The 
resistivities are the reciprocals of MC averaged conduc- 
tances. The normalized magnetization (Af ) is obtained 
from the spin structure factor S{k), ai k — (0,0).— 



III. RESULTS AND DISCUSSION 
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To start the discussion of results, the original state 
without any anisotropic contribution is simulated as a 
reference. The resistivities along both the x and y di- 
rections [px and py) are calculated as a function of tem- 
perature (T), as shown in Fig. 1. As expected, px and 
Py are almost identical in the whole T range. The small 
differences between px and py are from statistical fluctu- 
ations during the MC simulation, and these differences 
should converge to zero with increasing MC simulation 
times. With this set of parameters, both px and py show 
a MIT with increasing temperature at Tmi '^ 0.045io, 
which is the same approximate location as our estima- 
tion for the Curie temperature (Tc), according to the 
M — T curve. For a typical manganite with a MIT un- 
der zero magnetic field, Iq is roughly estimated to be in 
the range 0.4 - 0.5 eVJ^i22 Thus, Tmi ^ 200 - 260 K in 
agreement with bulk measurements. Therefore, the set 
of parameters (Jaf — 0.09, A = 1.2) used here is suitable 
to describe typical manganites, such as Lai_2:Caa;Mn03. 

In the following, we will apply the aforementioned 
three anisotropic interactions one by one into the model 
simulation to clarify their respective roles. First, let us 



FIG. 2. (Color online) MC simulated resistivities (triangles) 
and magnetization (dots) for strained lattices, as a function of 
T. (a) Only the anisotropic BE interaction is here considered 
[At = 0.1, Aj = 0, and Aq = 0). (b) Only the anisotropic 
SE interaction is considered [At — 0, Aj = 0.1, and Aq — 0). 
To maintain the presence of a metal-insulator transition J^p 
must be slightly reduced to 0.086. (c) Only the strained JT 
distortion is considered [At = 0, Aj = 0, and Aq — 0.01). 



consider the anisotropic DE interaction. For this pur- 
pose, At is set to 0.1 while other parameters are kept 
the same as the original ones. In other words, the DE 
hopping amplitude along the y direction is made 10% 
larger than that along the x direction, because of the 
presence of more straight Mn-O-Mn bonds along the y 
axis. The resistivities and magnetization of this strained 
lattice are shown in Fig. 2(a) as a function of T. px 
shows a MIT similar to the original one while py is now 
considerably suppressed in magnitude. Thus, a high de- 
gree of anisotropic resistivities can be obtained using an 



anisotropy At in the hoppings which is only 0.1. Inter- 
estingly, although the difference between p^ and Py are 
substantial, the differences in the Tmi's shown in p^ and 
Py are not obvious in this At = 0.1 case. Comparing 
with the original one, Tq and Tyii actually simultane- 
ously raise to ^ 0.055io due to the increase in ig. 

Next, the anisotropic SE is taken into account. At is 
restored to 0, while Aj is set to 0.1. In this case, J^p 
has to be weakened slightly to preserve the presence of 
an MIT, otherwise the system becomes insulating in the 
whole T range if J^p remains at 0.09. Thus, for this 
case the new values J\^ = 0.085 and J^p = 0.0935 are 
adopted. The MC simulated resistivities and magnetiza- 
tion for this strained lattice are shown in Fig. 2(b), as a 
function of T . The Tmi remains isotropic and coincides 
with Tc ~ 0.05io- In contrast to the DE case, the dif- 
ferences between p^ and py are much smaller, especially 
below Tc (or Tmi): Py is only slightly lower than px above 
0.04io, and they are almost identical below 0.04to. Then 
we conclude that the effect of an anisotropic SE is much 
weaker than the case of an anisotropic DE, when their 
anisotropic ratios are the same. 

Finally, it is necessary to address the effect of 
anisotropics in the JT sector, for completeness. In a dis- 
torted oxygen octahedron, the two eg orbitals are not 
degenerate anymore. For instance, when the lattice con- 
stants along the x and y axes are different, as in the 
Ward et aVs strained manganites thin films, < Q2 > ^s 
no longer zero. This nonzero < Q2 > mode induces an 
orbital-state "preference" over the whole lattice. With 
At = 0, Aj = 0, and Aq = 0.01, the MC simulated re- 
sistivities and magnetization are shown in Fig. 2(c), as a 
function of T. Similar to the case of an anisotropic DE, 
there is now a substantial difference between p^ and py. 
In addition, the Tmis of the px and py curves becomes 
anisotropic: the lower resistivity curve has a higher Tmi, 
in agreement with the experiments4S. 

To further clarify the anisotropic resistivities observed 
here, the relative percentage difference (6) between px 
and Py (defined as 6 = {px — Py)/ Py x 100%) is calculated 
for each of the three cases discussed above, as shown in 
Fig. 3(a). For the original isotropic and the Aj = 0.1 
cases, the values of 5 are very small (< ±20%) in the 
whole temperature range, as expected from Figs. 1 and 
2(b). In contrast, for the At — 0.1 and Aq = 0.01 cases, 
the situation is different. With increasing T from low 
temperatures, 5s first increases. After each case reaches 
a robust peak of 200 — 300%, then they decrease with 
further increases in T. Interestingly, for both these two 
cases, the corresponding Ts of the peaks found in 5 are 
slightly lower than the corresponding Tcs and TmiS, in 
agreement with the experimental results.— 

To understand the physical mechanism leading to 
the anisotropic resistivities the orbital properties of the 
strained states, characterized by the average values of the 
pseudo-spin orbital operator < t^ >, are also calculated, 
as shown in Fig. 3(b). The occupation difference between 
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FIG. 3. (Color online) (a) Relative percentage differences 
between p^c and py, as a function of T, for the cases indicated. 
O stands for original isotropic, (b) The average values of 
< Tx > , as a function of T. (c) Sketch of the effect of strain on 
the orbitals. To better distinguish the orbital-leaves along the 
X and y directions, here < Tj; > / < n > is magnified to 0.1. In 
the sketch, the overlap of electronic clouds becomes stronger 
along the y direction and weaker along the x direction. Thus, 
the conductances, which are in proportion to the overlaps, 
become anisotropic. Note that the real overlaps are indirect 
and mediated by oxygens (not shown here). Inset: sketch of 
an oxygen octahedron's in-plane distortion. 



portion to < Ta; >. The values of < r^; > for the original 
isotropic case fluctuate around zero in the whole T range 
analyzed, implying that the weights of the \'iy^ — r^ > 



and |3x^ 



> orbitals are equal, as expected by sym- 



metry. For the Aj = 0.1 cases, < t.^ > still remains very 
small, implying that the anisotropic Jaf used is not rele- 
vant to affect substantially the orbital composition of the 
state. In fact, both these two cases give rise to (almost) 
isotropic resistivities. In clear contrast, for the At — 0.1 
and Aq = 0.01 cases, finite values for < r^, > are ob- 
served at high T, which are gradually suppressed by the 
FM transitions with decreasing T. For the Aq = 0.01 
case, the finite < Ta, > is mainly caused by the JT distor- 
tion, which remains finite at low T as long as the lattice 
is anisotropically distorted. However, the finite < Tx > 
for the At — 0.1 case is caused by the enhanced DE pro- 
cess along the y direction. Namely, it is a DE mediated 
polarization of the orbital occupancy. Thus, for the fully 
FM state at low T, this DE mediated orbital rearrange- 
ment is largely suppressed to near zero, which is different 
from the results obtained for the JT distortion case. In 
summary, in our simulation the large anisotropy of the 
resistivity emerges in those cases where there is an un- 
balanced in the orbital state population, as sketched in 
Fig. 3(c), although the value of 5 is not linearly depen- 
dent on < Tj; > in the whole T range. In simple terms, 
the orbitals that increase their overlaps due to strain are 
now more populated than the other ones. 

Note that all the above simulations were carried out on 
relatively small 8x8 clusters using clean- limit models and 
still the anisotropy observed is comparable to that found 
experimentally. This implies that clean-limit strained 
manganites can be as anisotropic as phase separated com- 
pounds. Therefore, the LPCMO phase separation and 
classical percolation at the (sub) micrometer scale does 
not appear to be essential to obtain highly anisotropic 
resistivities, but of course in the clean limit the strain 
induced by substrates must be sufficiently large to gen- 
erate a ^t = 0.1 as used here, while for phase-separated 
compounds this anisotropy arises from phase competi- 
tion. Thus, the high anisotropic resistivities should be a 
general properties of manganites and even other complex 
oxides, as long as the bond lengths/angles are tuned to 
be sufficiently anisotropic by strain. Further experimen- 
tal studies on strained oxide films are needed to verify 
our results. 

However, it is important to clarify that in the par- 
ticular case of large scale phase-separated manganites, 
the classic percolation mechanism can certainly also con- 
tribute to the anisotropic resistivities if the shapes of the 
FM metallic clusters become anisotropic, as suggested in 
Ref. 40. In fact, our model can also qualitatively explain 
the formation of anisotropic FM clusters. To study an 
individual phase-separated FM cluster embedded in the 
AFM COI matrix, the ground state energies of FM lat- 
tices with open boundary conditions can be calculated 
directly. For simplicity, all Q2i are set to be uniform 
(and equal to < Q2 >) and all spins are aligned to be 
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FIG. 4. (Color online) Energy differences of the ground states 
of Lx X Ly — 900 lattices, as a function of L^ ■ The energy of 
the Lx = 30 lattice is set as the reference point, (a) Results 
obtained with the anisotropic DE interactions, (b) Results ob- 
tained with the strained JT distortions, (c) Results obtained 
with both the anisotropic DE interactions and strained JT 
distortions simultaneously active. The values of {At, Aq) in 
(c) are simultaneously stepped the same as in (a) and (b), re- 
spectively, (d) Sketch of the formation of an anisotropic FM 
cluster, according to (a-c). 



perfectly FM. Then, the shape of the FM clusters can be 
determined by varying the lattice's shape but keeping a 
constant lattice area. For instance, the energies of lat- 
tices with the same area size {L^ x Ly = 900, L^ and 
Ly are side lengths along the x and y axes, respectively) 
are shown in Fig. 4(a-c), as a function of L^. The en- 
ergy of the 25 X 36 lattice, which is elongated along the 
y direction, can be obviously more stable than that of a 
30 X 30 one when At > 0.3, or Aq ^ 0.04, or At > 0.2 
and Aq ^ 0.02 simultaneously. This process is qualita- 
tively sketched in Fig. 4(d). FM clusters with other sizes 
(e.g. Lx X Ly = 576 and L^x Ly = 1764) have also been 



tested, reaching the same conclusion. Thus, it is reason- 
able to expect similar effects when FM clusters expand 
to the (sub)micrometer scale, although our microscopic 
model can not be directly used on such large lattices with 
the currently available computational capabilities. 

Finally, it is important to estimate how large should be 
the lattice mismatch required for the highly anisotropic 
resistivities observed here to appear in strained mangan- 
ites with nanoscale phase separation or in the case of 
a bicritical phase diagram. According to the well-known 
Harrison's formula,"*^ the DE hopping t and SE exchange 
Jaf can be estimated to be in proportional to r^^ and 
r~^^ (r is the Mn-0-Mn's bond length), respectively. 
Thus t/JAF is in proportion to r^. To obtain the val- 
ues At = 0.1 and Aj = 0.1 used in our simulations, the 
required lattice mismatch is about 1.4%. Similarly, by 



comparing experimental data (r; — rg ~ 0.6 A, where ri 
and Tg are long and short bonds, respectively)^ and the- 
oretical parameters (A|(52| — 1.5)— for the JT distortions 
in RMnOs, the parameter Aq — 0.01 used here is esti- 
mated as ~ 0.45%. It should be noted that the required 
strain (lattice mismatch in real films, or (At, Aq) in our 
simulations) depends on the particular materials under 
study (or, equivalently, the actual values of the param- 
eters (Jaf, A) in our simulations). The anisotropies are 
more sensitive to strain when the system moves closer 
to the phase boundary between the FM and COI phases. 
With this idea in mind, it is natural that the anisotropies 
of LPCMO films can be notorious even if the lattice mis- 
match is small in average (0.2 — 0.3% in Ward et aVs 
experiments), because LPCMO is precisely at the FM- 
COI phase boundary. According to our simulations, the 
highly anisotropic resistivities are also expected in other 
strained CMR manganites films (even without phase sep- 
aration), although the required strain might be somewhat 
larger than for the LPCMO case. 



ulation shows that the highly anisotropic resistivities 
are associated with an unbalanced in orbital popula- 
tions which is driven by the anisotropic double-exchange 
and anisotropic Jahn- Teller distortions. In contrast, the 
anisotropic superexchange was not found to be a domi- 
nant driving force for the anisotropic resistivities. The 
observed high anisotropic resistivities in our simulation 
did not rely on phase separation at the (sub) microscopic 
scale. Therefore, it is expected that this anisotropic state 
could be realized in a variety of manganites and other 
complex oxides as well, if a sufficiently large lattice mis- 
match can be achieved in the growth of the mangan- 
ite films. In addition, for the particular case of phase- 
separated manganites, our model investigations suggest 
that the anisotropic double-exchange and strained Jahn- 
Teller distortions could indeed reshape the ferromagnetic 
clusters, thus inducing an anisotropic percolation and 
concomitant anisotropic resistivity that further enhances 
these effects. 
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IV. CONCLUSIONS 

In conclusion, the high anisotropic resistivities of 
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